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Abstract. The density response of a doped Mott-Hubbard insulator is discussed starting from 
the t-J model in a slave boson \/N representation. In leading order 0(1) the density fluctuation 
spectra N(q, co) are determined by an undamped collective mode at large momentum transfer, in 
striking disagreement with results obtained by exact diagonalization, which reveal a very broad 
dispersive peak, reminescent of strong spin-charge coupling. The 1 /N corrections introduce the 
polaron character of the bosonic holes moving in a uniform RVB background. The resulting N(q, co) 
captures all features observed in diagonalization studies, fulfills the appropriate sum rules, and apart 
from the broadening of the collective mode shows a new low energy feature at the energy %J + 8t 
related to the polaron motion in the spinon background. It is further shown that the low energy 
structure, which is particularly pronounced in (tz,Q) direction, describes the strong renormalization 
and anomalous damping of the highest bond-stretching phonons in La2_ A Sr l Cu04. 

Dedicated to Ferdinando Mancini on the occasion of his 60th birthday 
published in AIP Conference Proceedings 695: "Highlights in Condensed Matter Physics", 
ed. A. Avella et al. (AIP, New York, 2003) 

INTRODUCTION 

High-temperature superconductors are doped Mott-Hubbard insulators, therefore the 
density response is expected to be very different from that of weakly correlated Fermi 
systems. The low-energy density response, i.e. in the frequency range inside the Mott- 
Hubbard or charge transfer gap, is proportional to the doping, i.e., described by transi- 
tions within the lower Hubbard band. For a complete understanding of the physics ruling 
these systems both the spin and the charge response must be considered in the full mo- 
mentum q and frequency CO domain. The spin response has been and still is extensively 
explored, mainly stimulated by a wealth of experimental data provided by neutron scat- 
tering. A similarly powerful experimental tool is missing for the charge response, with 
the exception of optical conductivity, which however provides information only in the 
limit q — * 0. In the absence of solid experimental information about the detailed struc- 
ture of N(q, co), the best we can do to test a theory is the computer experiment. This 
provides for small clusters an unbiased answer how Af(q, co) looks like, e.g., for the t-J 
modelling, and how this quantity depends on the key parameters, namely the magnetic 
exchange interaction J and the doping 8. 

In one dimension the Hubbard physics is basically simple (although not mathemati- 
cally), as it is characterized by charge and spin separation, which implies that the density 
response is basically that of non-interacting spinless fermions, i.e., showing vanishing 
excitation energy at q = and 4k f. This is even true for the ID t-J model away from 



7 = and 2t, where the model is not exactly solvable! 1]. As we shall see the 2D model 
relevant for the cuprates shows a very different behavior compared to the ID case. 

Exact diagonalization studies for the t-J modelyj, |2J] have revealed that the dynam- 
ical density response N(q, ft>) for the 2D model is characterized by strong spin-charge 
coupling. These calculations show several features unexpected from the point of view of 
weakly correlated fermion systems which have to be explained by theory: (i) the strong 
suppression of low energy 2kp scattering in the density response, (ii) a broad incoherent 
peak at high energy (several t) whose shape is rather insensitve to changes of hole con- 
centration and exchange interaction J, (iii) the very different form of N(q, ft)) compared 
to the spin response function 5(q, ft)), which share common features in weakly correlated 
fermionic systems, and finally (iv) there are at low energy (~ J) some features that do 
depend on doping and J. It is worthwhile to note that finite temperature diagonalization 
studiesO] show only weak temperature dependence for T < 0.3t even at low energy. 

While considerable analytical work has been done to explain the spin response of 
the t-J model only few authors analysed N(q, ft)). Wang et al.[4] studied collective 
excitations in the density channel and found sharp peaks at large momenta corresponding 
to free bosons. Similar results were obtained by Gehlhoff and Zeyher[5] using the X- 
operator formalism and by Foussats and Greco||a] using a path integral representation 
for X-operators. Lee et obtained a broad incoherent density fluctuation spectrum 
by considering a model of bosons coupled to a quasistatic disordered gauge field, while 
Jackeli and Plakida[8] employed the memory function approach. 
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FIGURE 1. Comparison of N(q, <x>) obtained by the slave-boson theory 1 9] (solid lines) with diagonal- 
ization data| 1] for a periodic 4x4 cluster with J/t = 0.4 and doping 8 = 0.25. The dashed line in the 
(n, n) spectrum indicates the 5-function collective peak obtained when polaron effects are neglected. 

Starting from a slave-boson representation we show that the essential features ob- 
served in the numerical studies can be obtained in the framework of the Fermi-liquid 
phase of the t-J model at zero temperature Ja]. Our main findings are: (i) at low mo- 
menta the main effect of strong correlations is to transfer spectral weight from particle- 
hole excitations into a pronounced collective mode. Because of the strong damping of 




this mode (linear in q) due to the coupling to the spinon particle-hole continuum, this 
collective excitation is qualitatively different from a sound mode, (ii) At large momenta 
we find a strict similarity of N(q, ft)) with the spectral function of a single hole moving 
in a uniform RVB spinon background. In this regime iV(q, ft>) consists of a broad peak at 
high energy whose origin is the fast, incoherent motion of bare holes, (iii) The polaronic 
nature of dressed holes leads to the formation of a second peak at lower energy, which 
is more pronounced in (tt,0) direction in agreement with diagonalization studies [jl 0]. 

We show that the anomalous renormalization of certain phonon modes as observed in 
inelastic neutron scattering provides a sensitive test of the pecularities of the low-energy 
density response. In this contribution we shall analyse the strong, doping-dependent 
renormalization of the highest breathing phonon modes which is a generic feature in the 
high-r c compounds. 



SLAVE BOSON THEORY OF DENSITY RESPONSE 

We start from the Af-component generalization of the slave-boson t-J Hamiltonian ifioi Ell 
H t j = H t +Hj, which is obtained by replacing the constrained electron creation operators 

H ' = -| I (f^hifja+h.c), (1) 
ly <i,j>a 

H J = i E ftafio'f^fja(l-hthi)(l-hpj), (2) 
<i,j>aa' 

where ff a is a fermionic (spinon) operator, o = 1, — ,iV" is the fermionic flavor index, 
and hj denotes the bosonic holes. These operators obey standard commutation rules, yet 
the number of these auxiliary particles must obey the constraint Y,afio fia+ hfhi = N/2. 
The original t-J model is recovered for N = 2. 

The slave boson parametrization provides a straightforward description of the strong 
suppression of density fluctuations of constrained electrons through the representation 
of the density response in terms of a dilute gas of bosons. A common treatment of model 
(1) is the density-phase representation ("radial" gauge llllh ) of the bosonic operator 
hi = r z exp(/0/) with the subsequent 1 /^-expansion around the Fermi-liquid saddle 
point. While this gauge is particularly useful to study the low energy and momentum 
properties, it is not very convenient for the study of the density response in the full co 
and q space. Formally the latter follows in the radial gauge from the fluctuations of rf. 
However, if one considers for example convolution type bubble diagrams, one realizes 
that their contribution to the static structure factor is correctly of order l/N, but is not 
proportional to the density of holes 8 as it should be. According to Arrigoni et al lfl2ll 
such unphysical results originate from a large negative pole contribution in the {r^r^) a 
Green's function of the real field r, which is hard to control by a perturbative treatment 
of phase fluctuations. 

We follow therefore Popov[13] using the density-phase treatment only for small 
momenta q < qo, while keeping the original particle-hole representation of the den- 



sity operator, h + h, at large momenta. More precisely /?,■ = r, exp(z'0 7 ) +&,-, where b{ = 
^| q i >go /z q exp(zqR i ). The cutoff qo is introduced dividing "slow" (collective) variables 
represented by r and from "fast" (single-particle) degrees of freedom bj and bf '. As 
discussed by Popov lfl3ll this "mixed" gauge is particularly useful for finite temperature 
studies to control infrared divergences. We start formally with the "mixed" gauge and 
keep only terms of order 8 and 1 /N in the bosonic self energies. In this approximation 
our zero temperature calculations become quite straightforward: The cutoff qo < 8 ac- 
tually does not enter in the results and we arrive finally at the Bogoliubov theory for a 
dilute gas of bosons moving in a fluctuating spinon background. 

The Lagrangian corresponding to the model (1) is then given by (the summation over 
o is implied) 

L = L{fta(^-^)fio+bt(^-il b )b^H t+ Hj 

{ftafio + in + bt)(r i + bi)-^), (3) 

H < = ~% L f&Mbpt + rtrj + rjbt + bfn) +h.c. (4) 

Here the X field is introduced to enforce the constraint, and \Xf,\l\, are fixed by the 
particle number equations (rif) = y(l — 8) and (rf + b^bi) = j8, respectively. The 
uniform mean field solution r ; = tq leads in the large N limit to the renormalized 

narrow fermionic spectrum ^ = —zt% — llf, with t = J% + t8, %, = j(cosk x + cosk y ), 
X = Y.a(fj^fja)/N, and z = 4 the number of nearest neighbors. In the N = °° limit 
Xoc — 2/n 2 is given by that of free fermions, while for the original t-J model its 
value should be larger full due to Gutzwiller projection. In the following % = \% x 
will be used. Distinct from the finite-temperature gauge-field theory of Nagaosa and 
LeeH the bond-order phase fluctuations acquire a characteristic energy scale in this 
approach[4], and the fermionic ("spinon") excitations can be identified with Fermi- 
liquid quasiparticles. The mean field spectrum of bosons is co q = 2zXt(l — y q ). Thus 
the effective mass of holes ml l/t is much smaller than that of the spinons. 

Due to the diluteness of the bosonic subsystem, 8 1, the density correlation func- 
tion Xq,co = (8n h 8n h )q a is mainly given by the condensate induced part which is repre- 
sented by the Green's function ((/>+ + b- q )(b q + /?l q ))co for q > qo, and 2(r_ q r q ) C() for 
q < qo, respectively: 

Xqco ^ y rl (((&+ + b- q ) (b q + bl q )) q>qo + 2{r- q r q ) q<q ^j . (5) 

The l/N self-energy corrections to these functions are calculated in a conventional 
wayyjlQlJ] expanding r,- = (ro\/N + (8r)i) /a/2 and considering Gaussian fluctuations 
around the mean field solution. Neglecting all terms of order 8/N and q^/N, only one 
relevant 1 /N contribution remains which corresponds to the dressing of the slave-boson 
Green's function by spinon particle-hole excitations. Within this approximation and at 



zero temperature no divergences occur at low momenta, thus one can take the limit 
qo — > 0. The final result for the dynamic stucture factor (normalized by the hole density) 

N^co = -Im^COqa+S^P -jUfr)/D q ,a>), (6) 

D q , a = (^+4!^ -/**) - (ofl-A^) 2 . (7) 

The origin of the contribution 

c(i) , 2/( 1 + n 2) 2 -n \ 

S y q , ! a = ztr ( — n 3 ) q(0 , (8) 

U m -ZtL-z ? (?k + ?k+q) , (9) 

k ''k+q _ w _ 60 ~ * u 

is the indirect interaction of bosons via the spinon band due to the hopping term (which 
gives in (4)) and due to the coupling to spinons via the constraint field A . The latter 
channel provides a repulsion between bosons, making (ft) = 0) positive and therefore 
ensuring the stability of the uniform mean-field solution. 

The l/N self energies S^ 1 /^) and A^ 1 /^ are essentially a single boson property. 
They are given by the symmetric and antisymmetric combinations (with respect to 
ft) + i0 + — > — ft) — /0 + ) of the self energy 

4,2° =15 E (^'-q) 2 ^ +k _ k ,(£0 + ^ k -^). (10) 

iv lkl<^<lk'l 

Here G q (ft)) = (ft) — ft> q — Eq 1 ^ + lib) 1 is the Green's function for a single boson 
moving in a uniform RVB background. Although in the context of l/N theory the 
G° function in (5) should be considered as a free propagator, we shall use here the 
selfconsistent polaron picture for a single hole[ 16]. This is crucial when comparing the 
theory for N = 2 with diagonalization studies. Finally, the constants a and in (3) are 
given by (1 — tr^/t) and S^' N ^(co = q = 0), respectively. The parameter rjj in EqEJ 
which formally corresponds to the condensate fraction in our theory, is determined 
selfconsistently from = 8 — Y^q^ohq. The momentum distribution hq = (b^bq) is 
calculated from the corresponding bosonic Green's function for finite hole-density. The 
interactions implied in the polaron formation lead to a significant reduction of the 
number of bosons in the condensate[9|], and affects the balance of the selfenergies 
and S^/AO i n Eq. (7). 

The structure of iV(q, ft)) (6) in the small ft),q limit is mainly controlled by the 
interaction of bosons represented by the term (°c Tq), while the internal polaron 
structure of the boson determined by S^ 1 /^) is less important. iV(q,ft>) consists of a 
weak spinon particle-hole continuum with cutoff vpq, and a very pronounced linear 
collective mode which nearly exhausts the sum rule. The velocity of this mode is always 
somewhat smaller than the spinon Fermi velocity, v s <vf— zt, which implies a strong 
damping ft) (or q) of this mode (Fig 13). 




FIGURE 2. Density fluctuation spectra N(q, to) for//? = 0.4 and 8 = 0.15 along (n,0) (left) and (tt, tt) 
directions (right). Energy in units of t. 
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FIGURE 3. Low energy density response at (k,0) for J = 0.3 and 8 = 0.18. The peak centered at 
~ (%J + 8t) due to the polaron motion appears in the l/N order on top of the spinon particle-hole 
continuum with edge at energy z{%J + 8t). Small oscillations in the data result from the discretization 
of the energy scale in the numerical solution and have no physical meaning. 




The density response N(q, co) (Fig. 2) at large momenta, q > 5, which we can compare 
with diagonalization results (Fig.l), is dominated by the properties of the single boson 
self energy S^/ N \ The calculated density response of the t-J model has three character- 
istic features on different energy scales: (i) The main spectral weight of the excitations at 
large momenta is located in an energy region of order of several t. This high energy peak 
is very broad and incoherent as a result of the strong coupling of bosons to low-energy 
spin excitations. The position of this peak and its shape are rather insensitive to the ratio 
J ft < 1 in agreement with conclusions of |2B - This is simply due to the fact that the 
high-energy properties of the t-J model are controlled by t. (ii) The theory predicts also 
a second peak at lower energy which is more pronounced in the direction (tt,0) (see 
also Fig. 3), while its weight is strongly suppressed for q near (tz, k). The origin of this 
excitation is due to the formation of a polaron-like band of dressed bosons. The relative 
weight of this contribution increases with 7 as a result of the increasing spinon band- 
width, (iii) In addition there is the spinon particle-hole continuum which is generated by 
with relatively small weight 8 2 ). At (tt,0) the high energy cutoff of the spinon 
continuum is at z(xJ + 5?) (Fig. 3), while the polaron peak is at about {%J '+ 8t). We 
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FIGURE 4. Comparison of static structure factor N(q) (solid line) calculated for J — 0.4 with the result 
for a Gutzwiller projected wave functionl 17| (squares) and spinless fermions (dash-dotted) for 8 — 0.213. 
The comparison with the result obtained for non-interacting fermions with spin (dashed line) shows the 
large reduction of the density response due to the constraint. 



note that the polaron peak in iV(q, Q)) at optimal doping is in the same energy range as 
the high energy longitudinal optical phonons in cuprates, and as we shall discuss below 
leads to a strong damping and anomalous renormalization of these modes for certain 
wave numbers. 

Figure 4 provides a comparison of the static charge structure factor N(q) with that of 
a Gutzwiller projected wave function and with that for uncorrelated electrons. The latter 
comparison shows the strong suppression of the density response due to correlations. 
Finally we note that our calculated structure N(q) factor satisfies the sum rule for 
constrained fermions = 5(1 — 5) within 1% [9]. 



Inelastic neutron scattering experiments on high-r c superconductors have shown that in 
particular the highest energy longitudinal optic phonon branch near (tt,0) softens and 
broadens strongly as holes are doped into the insulating parent compound. Whereas 
the corresponding breathing mode at (n, n) shows much smaller softening and no 
anomalous broadening (Fig O . This effect appears to be generic for cuprates and de- 
tailed neutron scattering studies have been reported for La2 x Sr x Cu04 111 81 
YBa9Cu^Oft +y ||19Ll22Ll23ll and also for Bi-based cuprates|24]. This certainly shows that 
there are phonons which couple rather strongly to the charge carriers that contribute 
to superconductivity, and one may ask what is special about the interaction of strongly 
correlated electrons and phonons. 

The renormalization of these phonons can be calculated in the framework of the t- 
J model, since the breathing oxygen motion of these modes (see FigO modulates the 
energy of a hole in a Zhang-Rice singlet state, and therefore couples directly to the 

t 1 

density of doped holes. Expanding the Zhang-Rice energy Ezr = 8;f§ with respect to the 
oxygen displacements u l a = u a (Rj + 8®), oc =x,y, of the four O-neighbors at Ri + 8° 



RENORMALIZATION OF BREATHING PHONONS 
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FIGURE 5. Doping dependence of the high energy bond-stretching phonons along (£,0,0) and 
(£, £,0) directions for undoped La2CuC>4 (open circles) and Lai.gsSro.isCuCU (filled circles) as obtained 
by neutron scattering. The lower panel gives the line width at half maximum. The displacements of oxygen 
ions (a) for the q = breathing phonon distortion and (b) for the (it,0) half -breathing mode are 

indicated by arrows in the right panel. Here (n, 0) corresponds to (0.5,0,0) reduced lattice units, (reprinted 
from Pintschovius und Reichardt lfl9ll with permission from World Scientific) 
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FIGURE 6. Calculated phonon spectral function B p i, for (n,Q) and (n,n) breathing phonons for 
8 = 0.15, t = 0.4 eV, / = 0.12 eV and £, = 0.25 (solid lines) compared to undoped system (dash-dotted 
lines) and the inelastic neutron scattering data|20| for LaigsSro 15G1O4. (reprinted from Khaliullin and 
Horsch[25] with permission from Elsevier) 



around the Cu-hole at 7?, yields the linear electron-phonon coupling ||2. 



-ph 



(11) 



—7 12 

We assume that the resonance integral obeys the Harrison relation t pc t «r , where 
ro is the Cu-0 distance, and obtain g = IEzr/^tq, i.e., g w 2eV/A. The lattice part of 
the Hamiltonian is determined by the force constant K m 25eV/A 2 for the longitudinal 
O-motion. Due to the structure of H e - p ) t the breathing modes couple directly to Xq.a>- 



We have studied the renormalization of the phonon Green's functions along (it, 0) and 
(it, it) directions 
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where (»q o is the bare phonon frequency, i.e. measured in the undoped parent compound, 
and a q = ^(sin 2 g x /2 + sin 2 q y /2). Based on the parameters of the pd-model we 
estimate for the dimensionless coupling constant £ = g 2 /tK ~ 0.3 — 0.5. The coupling 
constant a 2 vanishes at the T point and becomes maximal at the zone edges. This 
explains the marginal changes of renormalized phonon frequency (0 q at q = (0, 0) . The 
particularly strong increase of the renormalization for q along (it, 0) is a combined effect 
of the strong increase of a 2 and the low energy polaron structure in the density response. 
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FIGURE 7. Doping dependence of low-energy density response at (?r,0) (solid lines). As a conse- 
quence of the scaling of the polaron structure « (%J + St) there is a strong change in the renormalization 
and damping of the (tz,Q) optical phonon (dashed lines), which is at (Oq = Q.2t in the undoped system. 

Figure calculated for typical parameters for cuprate superconductors, shows the 
strong renormalization of the (tt,0) half-breathing mode for Lai .ssSro.isCuCU with a 
twice as large shift as for the (n, it) breathing phonon. The large damping of the (tt,0) 
phonon results from the hybridization with the large polaron peak in iV(q, (O) at this 
momentum and is consistent with the experimental data[20]. The phonon energies of 
the undoped parent compound (0 q o = 80(90) meV for (lt,0) and (it, it), respectively, 
are taken from Ref. il 1911 . 

The strong doping dependence of this effect, shown in Fig|7j is due to the scaling 
oc (%J + 8t) of the polaron peak position in N(q, 00). In particular our results imply 
that the large damping of the (lt,0) optical phonon should disappear at larger doping 
concentrations. 



SUMMARY 



We have outlined a 1/N slave-boson theory for the density response of the t-J model, 
which explains the data obtained by exact diagonalization. We demonstrated that the 
predicted low energy polaron structure in the density response, which is particularly 
pronounced along (tt,0), explains the anomalous doping induced line width and shift 
of the longitudinal planar (tt,0) phonon. The energy of the polaron peak is determined 
by the spinon energy scale, therefore we predict a nontrivial doping dependence for 
the phonon renormalization. In that respect further neutron scattering studies of the 
doping dependence of phonons would provide a sensitive test for the low energy density 
response as well as for the spin structure in the different doping regimes. 
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